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Abstract 

Empirical analysis of the foreign exchange market is conducted based on methods to 
quantify similarities among multi-dimensional time series with spectral distances in- 
troduced in [A.-H. Sato, Physica A, 382 (2007) 258-270]. As a result it is found that 
the similarities among currency pairs fluctuate with the rotation of the earth, and 
that the similarities among best quotation rates are associated with those among 
quotation frequencies. Furthermore it is shown that the Jensen- Shannon spectral 
divergence is proportional to a mean of the Kullback-Leibler spectral distance both 
empirically and numerically. It is confirmed that these spectral distances are con- 
nected with distributions for behavioral parameters of the market participants from 
numerical simulation. This concludes that spectral distances of representative quan- 
tities of financial markets are related into diversification of behavioral parameters 
of the market participants. 
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1 Introduction 



Developing and spreading of Information and Communication Technology (ICT) 
computerize our society all over the world. One can access information more easily 
than before 1980s, while one sometimes faces to occasions when one has to handle 
massive information beyond human information processing capacity. 

In early 1990s financial markets started to introduce electronic trading systems 
due to development of ICT. This brings that information generated in the finan- 
cial markets almost becomes detectable. It is even becoming possible to conduct 
algorithmic trading or automated trading by using real-time financial data. 
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The foreign exchange market is the largest financial market all over the world. 
Recently an increasing number of financial institutions are embracing the op- 
portunities and benefits in automated and algorithmic trading systems. These 
systems have been successively studied in order to avert risks due to explosion of 
information and several financial institutions even experimentally or commercially 
utilize them [1,2,3]. It is predicted that more becoming common more trading can 
be quickly conducted. Moreover they have an advantage to human traders from a 
viewpoint of fairness, accuracy, speed, and cost-effectiveness. If they excess criti- 
cal mass then human traders may not understand what happens in the financial 
markets. Therefore it is expected to become an important work to construct a sys- 
tem to quantify market states (states of market participants in financial markets) 
automatically since monitoring of financial markets of which speed is predicted 
to be beyond human information processing capacity is needed. From physical 
point of view practical methods should be considered. In order to construct such 
a system constructive utilization of high frequency financial data should be so- 
phisticated. 

Generally speaking massive data contain less information since one cannot un- 
derstand and handle them without computers. In order to understand them some 
kinds of information extraction techniques are proposed [4]. By application of the 
random matrix theory to correlation matrices of multi-dimensional time series 
some studies showed that market states reveal in both eigenvectors and eigen- 
values [5,6]. It is known that the bone structure consisting of larger eigenvectors 
is related to clusters of financial markets. This approach is based on methods to 
extract network structures of correlation network [7]. In fact correlations among 
financial commodities are related to market structure but multiple aspects to 
measure states of financial markets should be discussed. 

In this article we consider alternative method to quantify market states with 
spectral methods. Moreover the relation between spectral distance and behavioral 
parameters of agents is discussed through an agent-based model. Recent works 
show that agent-based models of financial markets are useful to understand price 
formation mechanism in the financial markets [8,9]. 

The rest of this article is constructed as follows. In Section 2 definition of spec- 
tral distance for multi-dimensional time series is presented. In Section 3 quotation 
rates and quotation frequencies are introduced as representative quantities of fi- 
nancial markets. In Section 4 the Jensen-Shannon spectral divergence and the 
Kullback-Leibler spectral distance are calculated by using high-frequency finan- 
cial data of the foreign exchange market. In Section 5 in order to check adequate- 
ness of proposal methods those spectral distances are calculated by an agent-based 
model. Section 6 is devoted into concluding remarks. 



2 Definition of spectral distances for multi-dimensional time series 

Consider M-dimensional time series Xj(k) (j = 1, . . . , M; k = 0, 1, . . . , L — 1) 
with a sampling period At. The power spectrum of Xj(k) is estimated as period- 



gram estimator with the Harming windowing function with width N (< L), 
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The periodgram estimator is calculated as 
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where the Nyquist frequency is / c = 1/(2 At). By using the normalized power 
spectrum, which is defined as 

Pj(fn,t) = N Ji - 77 — , (3) 

the spectral entropy is calculated as its Shannon entropy; 

N-l 

H(Pj, t) = - Pj(fn, t) logpj(f n , t). (4) 

n=l 

Note that a direct current component is ignored in this definition. Obviously when 
Xj(k) is a white noise it takes the maximum value since the power spectrum of it 
is a uniform function. Contrarily when Xj(k) is a sinusoidal function it takes the 
minimum value since the power spectrum of it is a Dirac's 5-function. 

In fact the spectral entropy indicates randomness of time series. High/low H 
means that Xj(k) is unpredictable/predictable. Suppose the situation to predict 
time series from historical observation. The mode for a normalized power spec- 
trum is the frequency where the normalized power spectrum has the maximum 
value: 



/_,■(*) = arg max{pj(f n ,t)}. (5) 
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If the spectral entropy is low then sinusoidal signal with fj(t) is an adequate 
model to predict the future time series. Contrarily if the spectral entropy is high 
then sinusoidal signal with fj(t) is unworkable to predict the future time series. 
Therefore the spectral entropy is related into predictability of the mode on a 
frequency domain. 

Actually according to the maximum entropy principle the best choice for the 
distribution Pj(f n ,t) is the one that maximizes Eq. (4) subject to the constraint 
the normalized condition 
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By using the Lagrange multiplier differentiating 
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with respect to Pj(f n ,t) and then equating the result to zero one obtains the 



uniform normalized power spectrum 

Pj(fn,t) = j^j. (8) 

Namely the best choice without knowledge is the uniform power spectrum, which 
means that time series are white noises. 

By using the spectral entropy the Jensen-Shannon spectral distance [10,11] for 
M-dimensional time series is defined as 
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where Y^jLi n j — 1 an d 7Tj > 0, one can quantify the similarity among a set of 
time series. The Jensen- Shannon spectral distance measures similarity of random- 
ness among time series. As an alternative spectral distance the Kullback-Leibler 
spectral distance defined as 

KL lm {t) = EWn,0logM^4 (10) 

n=l Pm\Jn,t) 

can be considered in analogy with relative entropy between normalized power 
spectra [12,13,11]. 



3 Representative quantities of financial markets 

What are representative quantities of financial markets? Since there is circu- 
lar causality between market participants and a broking system quotations are 
causes/results for broking systems/market participants, and transactions are re- 
sults/causes (see. Fig. 1). 
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Fig. 1. A pictorial illustration of circular causality of a financial market. 

For empirical investigation the database of the foreign exchange market in the 
CQG datafactory is used [14]. As representative quantities of the financial market 
quotation rates and quotation frequencies are detectable from this database. The 
data contain 20 currency pairs (AUD/JPY, AUD/USD, CAD/JPY, CHF/JPY, 
EUR/AUD, EUR/CHF, EUR/GBP, EUR/JPY, EUR/NOK, EUR/SEK, EUR/USD, 
NZD/USD, USD/CAD, USD/CHF, USD-/JPY, USD/NOK, USD/SEK, USD/ZAR, 



GBP/AUD, and GBP/JPY) which consists of 11 currencies; AUD (Australia Dol- 
lar), CAD(Canadian Dollar), CHF(Swiss Franc), EUR(Euro), GBP(United King- 
dom Pound), JPY( Japanese Yen), NOK(Norwegian Krone), NZD (New Zealand 
Dollar), SEK(Swedish Krona), USD(United States Dollar), and ZAR(South Africa 
Rand) . 

Quotation frequencies are defined as the number of quotations per unit time, 

Mk) = ^Cj(kAt; (k + l)At) (k = 0, 1, . . . ; j = 1, . . . , 20), (11) 

where Cj(ti, t 2 ) is the number of ask/bid quotations of the j-th currency pair from 
ti to t 2 , and At denotes the sampling period. At is set as 1 [min] throughout this 
analysis (f c = 0.5 [1/min]). The best ask/bid rates are defined as the minimum 
ask/maximum bid rates between AAt and (k + l)At. 
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where askjit) /bidjit) denotes ask/bid rates of the j-th currency pair at time t. 

In this article, on the basis of spectral methods quantifying dynamical structure 
of best quotation rates and quotation frequencies is conducted. 



4 Empirical analysis of the foreign exchange market 



As shown in Fig. 2 the Jensen-Shannon spectral distances for best ask rates and 
ask quotation rates fluctuate with rotation of the earth. The Jensen-Shannon 
spectral divergence for the ask quotation frequencies denoting JSa and that for 
the best ask rates denoting JSr are calculated as shown in Fig. 2. 
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Fig. 2. The Jensen-Shannon spectral divergence of 20 currency pairs for the ask quota- 
tion frequencies (a) and that for the best ask rates (b) for a period from 16th October 
2006 to 20th October 2006. 



In order to calculate a correlation between them a cross-correlation coefficient is 



introduced as 



Cra 



(JSrJS a )-(JSr)(JS a ) 



(13) 



/Var{JS R }^Var{JS A } 
It is estimated as C RA = 0.624966 for a period from 16th Oct 2006 to 20th Oct 
2006. As a result it is found that best ask rates and quotation frequencies are 
weekly and constantly coupled with each other. 

Furthermore it is confirmed that values of the Jensen-Shannon spectral distance 
are related to the Kullback-Leibler spectral distance matrices. A mean of the 
Kullback-Leibler spectral divergence, defined as 



1 



M M 
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is always greater than JS(t) ({KL)(t) > JS(t)). Furthermore it is found that 
JS(t) is proportional to (KL)(t) with a coefficient as 0.42 from the least square 
method for empirical results as shown in Fig. 3 (JS(t) oc (KL)(t)). 
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Fig. 3. The relation between a mean of the Kullback-Leibler spectral distance and 
the Jensen-Shannon spectral distance for 20 currency pairs. A straight line represents 
a slope (0.42) obtained from the least square method, and circles are obtained from 
empirical data. 



5 Agent-based model 



By introducing an agent-based model of a financial market in which N market 
participants exchange M financial commodities [11]. 

The agents have threshold dynamics like the Granovetter model [15]. Three invest- 
ment attitudes are codes as 1 (buying), (waiting), and -1 (selling). If the inner 
states of the i-th agent for the j-th financial commodity at time t (t = 0,1,2...) 
excesses buying threshold 0fj{t) or below selling threshold 6fj(t) then the i-th 
agent selects the investment attitude for the j-th financial commodity at time t 



i (<M0 >*§(*)) 

o (eg(f) < s y (f) < e§(t)) , 
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where ®ij(t) is assumed to be described within diversification of interpretation 
£i(t) and self-feeling a,ij(t): 

Qijit) =a ij (t)(x i (t)+Z i (tf), (16) 

where Xi(t) denotes perception of the i-th market participants. For simplicity 
the perception is assumed to be described as both exogenous and endogenous 
information; 

M I T 

Xi(t) = £ U\emi \<%(t)\)7r E r k (t - rAt) + 8i (t), (17) 

k=l 1 t=1 

where rj(t) = ^/NY^Li Uij(t) denote log returns, Si(t) exogenous information, 
Cik(-, •) functions to show attentions, and T(> 1) a time span to calculate a 
moving average. The attentions are assumed as Cik{x,y) = l/(x 2 + y 2 ). Rates 
Rj(t) and quotation frequencies Aj(t) of the j-th financial commodity at time t 
are respectively defined as 

/ \ 1 N 

R j (t + At) = R j (t)exp(r j (t)), Aj(t) = ^EW*)I ( 18 ) 



Relationships between spectral distance and behavioral parameters of the mar- 
ket participants are investigated. It is found that diversification of parameters 
of agents is related to the Jensen-Shannon spectral distance for behavioral fre- 
quencies as shown in Fig. 4 (a). Furthermore it is confirmed that a mean of the 
Kullback-Leibler spectral distance matrices is proportional to the Jensen-Shannon 
spectral divergence from numerical simulation as shown in Fig. 4 (b). 



CO 




Fig. 4. (a) The relationship between parameters and JS of behavioral frequencies. The 
diversification of is measured as an entropy of parameter distribution. Parameters 
are fixed as N = 2000, M = 20 and T = 1. 0fj(t), 0fj(t), and a^t) are sampled 
from uniform distribution. The ranges of thresholds are fixed as [—0.02, —0.01] and 
[0.01,0.02], and the range of aij{t) [01,02] is measured by entropy H a = log(a2 — ai); 
(b) The relation between a mean of the Kullback-Leibler spectral distance and the 
Jensen- Shannon spectral distance. A straight line represents a slope (0.42) obtained 
from the least square method, and circles are obtained from numerical simulation. 
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6 Conclusion 



Spectral distances can be adopted to quantity to measure the market states to- 
tally. As a result it is found that the spectral distances fluctuate with rotation 
of the earth and that best quotation rates are weakly coupled with quotation 
frequencies. The Jensen-Shannon spectral distance is proportional to a mean of 
the Kullback-Leibler spectral distance. 

It is expected that employment of both spectral methods and high-frequency fi- 
nancial data provides us deepen insights of financial markets with high resolution. 
Information extraction from massive data will become important to understand 
states of financial markets. 

This work was partially supported by a Grant-in- Aid for Scientific Research 
17760067) from the Japanese Ministry of Education, Culture, Sports, Science and 
Technology. 
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